PROCEEDINGS 

OF SCIENCE 



D semi-leptonic decay form factors with HISQ charm 
and light quarks 



H. Na™, C. T. H. Davies 6 , E. Follana c , P. Lepage^ and J. Shigemitsu 

"Department of Physics, The Ohio State University, Columbus, Ohio, USA 
Department of Physics and Astronomy, University of Glasgow, Glasgow, UK 
' Departamento de Fisica Tedrica, Universidad de Zaragoza, Zaragoza, Spain 
d LEPP, Cornell University, Ithaca, NY, USA 



E-mail: heena@mps . ohio-state . edu 



We present a program to study D semi-leptonic decay form factors with HISQ charm and light 
quarks. In this exploratory work, we study D s to 7J S , /V/ semi-leptonic decay on MILC coarse 
(20 3 x 64) gauge configurations with 2+1 dynamical flavors. We develop a new method to calcu- 
late /+(0) and fo{q 2 ) from scalar current matrix elements, < 7r|sc|Z) >, which does not require 
any operator renormalization. We also determine f+(q 2 ) and fo(q 2 ) from vector currents with a 
fully non-perturbative renormalization method. 



The XXVII International Symposium on Lattice Field Theory - IAT2009 
July 26-31 2009 

Peking University, Beijing, China 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa. it/ 



D semi-leptonic decay with HISQ action 



H. Na 



1. Introduction 

Understanding the D —> % and D — > K semi-leptonic decays is an important topic in Charm 
physics. Furthermore, lattice studies of D semi-leptonic decays are of particular interest for several 
other reasons as well. First of all, one can test lattice formulations of generic semi-leptonic decays 
by studying the D semi-leptonic case. For example, \V„b\ and \ V c b\ can be determined from B semi- 
leptonic decays, and these are very important parameters to understand CKM physics and test the 
Standard Model. However, simulating B mesons is, in general, more difficult than D mesons. In 
B semi-leptonic decays, the maximum momentum transfer, q 2 ^, is much larger than in D semi- 
leptonic decays. In fact, one can cover the entire q 2 region of D semi-leptonic decays on the lattice, 
while, for B semi-leptonic decays, that would be challenging. In addition, the charm sector has 
been studied extensively in experiments as well. Thus, we would like to test the lattice formula- 
tion in the charm sector. Moreover, studying D semi-leptonic decays can provide an independent 
determination of \V CS \ and \V c d\- 

The D semi-leptonic decay form factors, f+(q 2 ) and fo(q 2 ), can be defined as 

(ttiv^iz^) =r+(^ 2 )[^+^ -^^^^i+y-oC^ 2 )^^^^, (i.i) 

where the momentum transfer q^ = p% — p%- Note that 

/ + (<7 2 = 0)=/ (<7 2 = 0), (1.2) 

because of kinematics. This is the conventional definition of form factors, which is normally used 
in experiments. For lattice calculations, it is more convenient to define the form factors as 



(n\V»\D) = v^mD[v^yj|(£»)+^/L (£»)], (1.3) 

where v M = and p± = p% — {j>% ■ v)v^ . This definition is particularly useful for lattice calcula- 
tions in the D rest-frame (v = (1,0,0,0) and p ± = (0,p n )), because each form factor is determined 
from the temporal current matrix element or the spatial current matrix element respectively; 

^ E « )={ -^ [GeVl/2] > f ^ )=i ^vJ GeV " /2] - (L4) 

So, the traditional strategy is to calculate /j| and f± on the lattice, take the continuum and chiral 
limit, and then convert to /+ and /o using the following relations; 



Uq 2 ) = "f^ 2 [(m D - m x )f\\ (E a ) + (E 2 n - m 2 n )f^E n % (1.5) 
mo- m 2 

f + (q 2 ) = -^=lf^ En ) + {m D -E n )f L {E n )\. (1.6) 

We can, then, compare the lattice calculation to experiments. In this talk, we present a new strategy 
that allows us to calculate the relevant form factors from the scalar current (S) rather than the vector 
current (V M ) (Sec. 2 and 3), and the more traditional method using (V^) (Sec. 4). 
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2. Formalism 

We apply the HISQ action [[lj] for all valence quarks; light, strange, and charm quarks. As we 
use the HISQ action, there are several noticeable advantages. We can simulate relativistic charm 
quarks, since the HISQ action has small enough discretization errors [|J]. In addition, we can 
formulate a fully non-perturbative renormalization. 

Consider the following relation between the vector and scalar currents in the continuum; 

q»(V™ nti ) = (m c -m q )(S conti ), (2.1) 

where (Vn 0nti ) = (7t\V™ ntl \D), (s contl ) = (n\sc\D), and q represents a strange or light quark. The 
RHS of Eq. ^T] does not require any operator matching factor (the combination is renormalization 
group invariant). Thus, one can write down the same relation on the lattice with a operator matching 
factor Z just on the LHS; 

q»(Vf t )Z=(m c -m q )(S lat ), (2.2) 

where S ,at is the local scalar current. This relation is only valid when we use the same action for 
all valence quarks, and the action has enough chiral symmetry. If one uses different types of action 
for different quark species, then the RHS of Eq. |2j] is no longer RG invariant. This is because the 
mass renormalization factors for charm and light quarks would be different, in general. 



In the D rest-frame, Eq. 2.2 can be re-written 



(m D -E x ){Vo)Z,+p x ■ (V)Z S = (m c -m q ){S), (2.3) 

where Z t and Z s are the Z factors for temporal and spatial vector current respectively. We now 
omit a superscript 'laf on the currents. Note that Z t and Z s can be different in general. Using this 
relation, one can extract Z factors fully non-perturbatively. We will explain this in more detail in 
Sec. 4. 

The scalar current matrix element can be written in terms of a single form factor fo\ 

= rn^-ml 2 
m r — m„ 



One can simply derive this relation from Eq. |ET] and 2.1. Therefore, we can obtain fo(q 2 ) from 
(S), 

Atf)= 0*-",>ffl , (2 . 5) 
m D-K 

The numerator of the RHS is the same RG invariant combination as in Eq. |2.1[ , and the denominator 
is a mass difference of the bound-states. Therefore, we can calculate fo(q 2 ) with no need for 
operator matching. This is the most interesting feature due to the fact that we use the HISQ action 
for all valence quarks. In addition, using this method, we can determine / + at q 2 = 0, because of 



Eq. 1.2; 



M0)=M0) = ^^\ . (2.6) 
nip — ra£ q 2 =0 

/+(0) is a very important quantity, since we can estimate |V CJ | and \V cc i\ by combining /+(0) and 
experimental inputs. 
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3. Test run for f (q 2 ) and /+(0) from (S) 

We have performed an exploratory analysis with MILC configurations with 2+1 dynamical 
flavors. The configurations that we use are 300 coarse lattices (20 3 x 62) with sea quark masses 
uocrnii = 0.01 and uoam s = 0.05 We calculate charm and strange quark propagators using 
the HISQ action with quark masses am s = 0.0546 and am, = 0.66. For each configuration, we 
calculate propagators with a point source at two different time slices. We insert four discrete 
spatial momenta; p = (0,0,0), (1,0,0), (1, 1,0), and (1,1,1). The non-perturbatively determined 
renormalization factor e of the Naik term of the charm quark is —0.21 [^]. 

In this test run, we calculate fo(q 2 ) of D s — ► r\ s , Ivi semi-leptonic decay using Eq. [2~3[ We 
use the local scalar current, and the Goldstone taste for the D s and tj v . In the analysis, we perform 
Bayesian fits to the three-point function and the two-point functions of D s and r] s simultaneously. 
We use fit model functions for two and three-point functions, 

nexp nexpo 

C*(f) = 1 £af(e- E " t + e- E;! ^)+ £ (-l)^f (^Ve^^), 

nexp nexpo 

C?(0 = £ of (e-^' + e-Wr-t)) + £ (-1)^^' + (3-D 

nexp nexpo 

C 3 (t) = £ £ [^ e -* e -^) + (-l) r ^-^ e -W-') 

H=l 0=1 

where the states with a bar are the negative parity states, and T is the distance between the sources 
of the light and heavy mesons. We choose T = 19 for this test run. We can impose stronger 
constraints on the meson masses and amplitudes by taking larger fit time domains for the meson 
two-point functions. Note that the correlator of T] s at zero momentum consists of the positive parity 
states only. 

The result is shown in Fig. |[ The four data points in the figure represent p = (1, 1,1), (1,1,0), 
(1,0,0), and (0,0,0) from the left. We take an average over equivalent momenta to increase statis- 
tics. The result of p = (1,1,1) is accidentally very close to q 2 = 0, so in this particular test case, 
we do not need to worry about extrapolations to q 2 = 0. We expect to achieve much smaller errors 
in the future, since we can improve our calculation using random-wall source techniques, more 
statistics, and multiple T values. 



4. Test run for the form factors from (V^) with a fully non-perturbative matching 

We showed a new method that will lead to the most accurate determination of f + (q 2 = 0), 
however we are still interested in the calculation of the form factors using the traditional method to 
test the new method and estimate the q 2 dependence of /+. The traditional method was described 
in Sec. 1. 

The vector current and pseudoscalar operator in the continuum are 

V l ° cal = %(x)7nVM <° ca ' = %(x)75¥M, (4.1) 
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f from <S> 



§ 
2.4% 



1 2 

q 2 (QeV 2 ) 

Figure 1: fo(q 2 ) from (S). The errors are statistical errors only. The two numbers below the data points are 
the relative errors. 

which are local operators. However, we cannot simply use these operators on the lattice with 
the Goldstone taste mesons, because we use the staggered fermion formalism. In order to ob- 
tain non-zero matrix amplitudes, we need to consider non-local vector current or non-Goldstone 
pseudoscalar operators, for example 

rone—link 



v, 



ii 

mon—Gold 
% 



W q {x)y ll ^ q2 {x + ii) -► Vqiy)^ 8>i]V fc (y) 
Wq(x)ybj 5 \\r q2 (x) -► y q (y)[Yo75 ® 75%]*^ (y), 



(4.2) 



where *P is the fermion field in the spin-taste basis. In this study, we choose the spatial and temporal 
currents differently; 

(yone-iink} = (Goldstone T],|V^" n *|Goldstone 



(Vff 



local \ 



(Goldstone T] s \V Q local |non - Goldstone D s ) 



(4.3) 



We tested the temporal matrix elements with all possible combinations of Goldstone and non- 



Goldstone D s and r\ s , V Q one lmk , and V local . We obtained consistent results for all combinations, 



and the choice of Eq. |13| gives the best results. These combinations give the smallest statistical 
errors, because the taste splitting for the Ds is smaller than that of the r\ s . 

4.1 fo(q 2 ) 



In the traditional method, we need to calculate first fu and f± separately as explained in Eq. |1.4 . 
The results for fu and f± before the renormalization using Eq. are shown in Fig. ^ (a). 

After we get fo using Eq. T5 (green diamonds in Fig. |2| (b)), we need to estimate the operator 



matching factors. Taking into account that Z t and Z v can be different, we re-write Eq. |L5| with Z 
factors; 



Mq 2 



mi 



mi 



[{m D - m K )Z t f {l (E n ) + (El - m 2 n )Z s f x {E K )]. 



(4.4) 



In fact, we do not need to consider Z s for fo, since (E% — m|) is small relative to (mp — m n ) and Z s 



is close to one. Using Eq. [1J at p = (0,0,0), we can obtain Z t 

(5)1^(0,0,0) 



m, 



m K 



m D -E n (V<°™') 1^(0,0,0) 



1.037(5). 



(4.5) 
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Figure 2: /j| and /j_ from the traditional method are shown in (a). fo(q 2 ) without Z r and with Z r are shown 
in (b). fo(q 2 ) from (S) is also shown for comparison. Form factors are calculated from the same q 2 values, 
but data points are shifted in the figures for clarity. 



Using this Z t and Z s = 1, we get fo(q 2 ) (red squares in Fig. |2| (b)), which is consistent with fo(q 2 ) 
from (S). 

4.2 f + (q 2 ) 

We can also re-write Eq. L6 with Z factors for f+(q 2 ), 

f + (q 2 ) = -^=[Z t f\\ (E n ) + (m D -E n )Z s f ± {E n )}. (4.6) 

Z s now plays an important role for /+. When we fix Z t at p = (0,0,0), then we can obtain Z s by 

_ (m c - m s )(S) - {m D -E n )Z t (p = (0,0,0))(V Q to ^) 

p n . (y one-link^ 

since Z factors are not supposed to depend on momentum. The result is shown in Fig. || (a). As 
one sees, the errors of Z s are much larger than the error of Z t . This is because Z s is estimated from 
subtracting two positive quantities. As a result, the error of Z s is large, even though the relative 
error of each term in Z s is comparable to the relative error of Z t . The result of f+(q 2 ) is shown 
in Fig. | (b). This is a very preliminary result. We hope that we can address the best strategy to 
calculate f+(q ) with correct Z factors in the future. 



5. Summary and future plan 

We presented a new method to calculate form factors of D semi-leptonic decays. We showed 
how we calculate f+(q 2 = 0) and fo(q 2 ) with no need for operator matching in Sec. 2 and 3. This is 
a general approach applicable to any semi-leptonic decays. For instance, one can apply this method 
to K — > % semi-leptonic decay. If we have relativistic bottom quarks, we can calculate B — > % and 
B — > D*(D) semi-leptonic decay form factors as well, although one will still have the problem of 
large q 2 max . The important thing is that the action should have enough chiral symmetry and all 
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Figure 3: Z factors (a), and f+(q 2 ) from the traditional method (b). Form factors are calculated from the 
same q 2 values, but data points are shifted in the figures for clarity. 

valence quarks should be calculated from the same action. In this exploratory work, we showed 
that the HISQ action is quite promising for D semi-leptonic decay with this new method. It should 
be possible to obtain / + (q 2 = 0) with significantly smaller errors than the current published theory 
errors of ~10% [Q]. It is also possible that we can replace the most accurate determination of 
\V c d\ from neutrino and antineutrino scattering We plan to improve our calculation by applying 
random-wall source techniques, adding more statistics, using smaller lattice spacings, and working 
with multiple T values. We will also investigate the continuum and chiral extrapolations, and 
extrapolations to q 2 = to complete this project. 

In Sec. 4, we also presented the traditional method based on (V^) with fully non-perturbative 
renormalization. We confirm that this traditional method and the new method are consistent with 
each other. It gives larger errors, however it is still interesting. Using this method, one can get the 
q 2 dependence of /+. We are still exploring other strategies for non-perturbative renormalization 
of Vfi, which will reduce our current error on Z s . 
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